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• Background and Aims Plant growth and respiration still has unresolved issues, examined here using a model. 
The aims of this work are to compare the model's predictions with McCree's observation-based respiration 
equation which led to the 'growth respiration/maintenance respiration paradigm' (GMRP) - this is required to 
give the model credibility; to clarify the nature of maintenance respiration (MR) using a model which 
does not represent MR explicitly; and to examine algebraic and numerical predictions for the respiration: 
photosynthesis ratio. 

• Methods A two-state variable growth model is constructed, with structure and substrate, applicable on plant to 
ecosystem scales. Four processes are represented: photosynthesis, growth with growth respiration (GR), senes- 
cence giving a flux towards litter, and a recycling of some of this flux. There are four significant parameters: 
growth efficiency, rate constants for substrate utilization and structure senescence, and fraction of structure 
returned to the substrate pool. 

• Key Results The model can simulate McCree's data on respiration, providing an alternative interpretation to the 
GMRR The model's parameters are related to parameters used in this paradigm. MR is defined and calculated in 
terms of the model's parameters in two ways: first during exponential growth at zero growth rate; and secondly at 
equilibrium. The approaches concur. The equilibrium respiratiomphotosynthesis ratio has the value of 0-4, 
depending only on growth efficiency and recycling fraction. 

• Conclusions McCree's equation is an approximation that the model can describe; it is mistaken to interpret his 
second coefficient as a maintenance requirement. An MR rate is defined and extracted algebraically from the 
model. MR as a specific process is not required and may be replaced with an approach from which an MR 
rate emerges. The model suggests that the respiration:photosynthesis ratio is conservative because it depends 
on two parameters only whose values are likely to be similar across ecosystems. 

Key words: Respiration, growth, maintenance, model, senescence, recycling, photosynthesis, gross production, 
net production, respiration:photosynthesis ratio. 



INTRODUCTION 

There continues to be considerable discussion of respiration. 
This is focused partly on how to handle maintenance 
respiration (MR) when working within the 'growth 
respiration/maintenance respiration paradigm' (GMRP) 
(Amthor, 2000) and partly on the ratio (r) of respiration (R) to 
gross photosynthesis (P g ) in plant ecosystems, defined as 
r R P . There is concern as to whether r R P is constant or variable. 
The topic is important in addressing some climate change 
issues, such as short-term and long-term contributions of 
plant ecosystems to carbon (C) sequestration. On one hand are 
simplifying and perhaps optimistic holists, who would like to 
assume that r R:P is approximately constant over a range of 
species and conditions (e.g. Gifford, 2003). This might permit 
construction of simple ecosystem models applicable to 
climate problems. A recent example of this is Van Oijen et al. 
(2010), who, working within GMRP, with some assumptions, 
demonstrate that a more or less constant value of the ratio is 



to be expected, partly as a result of conservation of matter (stoi- 
chiometry). On the other are reductionists who might also con- 
sider themselves as realists. The latter, to which group I mostly 
adhere, make efforts to understand the detailed mechanisms of 
respiration which cause r R:P to fall below unity and, possibly, 
be rather variable (Cannell and Thornley, 2000; Thornley and 
Cannell, 2000). Amthor (2000) gives an excellent review of 
plant respiration, supplying a valuable historical perspective 
and also the current state of the art - there have been no substan- 
tive advances in the past decade. 

Mainteneance respiration is often considered to be a distinct 
portion of respiration, but it has long been regarded sceptically. 
Wohl and James (1942), remarked that 'It is a facile assumption 
that the energy change revealed by the continuous liberation of 
heat from mature tissues which are neither doing external work 
nor synthesizing appreciably, is a measure of the energy of main- 
tenance. The assumption will not stand closer inspection 

The construction and analysis of simple models can some- 
times be employed to illuminate important problems; this is 
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what is attempted below. A simple model without a mainten- 
ance component is described. Solutions are obtained, focusing 
on two important phases in the growth dynamics: exponential 
growth and the steady state. The first phase, exponential 
growth, is transient, but it often lasts sufficiently long to be 
characterizable, and it is amenable to experimentation. The 
second phase, the steady state, is important because many 
plant ecosystems may be approximately in a steady state, 
although it can be less accessible to experimentation. 

We define the specific MR rate (p m ; d _1 ) as being the 
specific respiration rate when the growth rate is zero, i.e. 
when the status quo is maintained. The growth rate can be 
zero in two cases: when the plant is in the exponential 
growth (eg) phase but is growing at zero rate; and when the 
plant is in the steady state (ss) (its mass has reached an 
asymptote). In both cases, the same expression is obtained 
for the specific MR rate, depending on three parameters of 
the model (Fig. 1). This permits a re-parameterization of the 
model, which can, with assumptions, be simplified to a form 
perhaps more acceptable to the holists. 

Current understanding of plant respiration was given 
impetus and direction by the findings of McCree (1970), and 
his famous equation: 



R = kP + cW 



(1) 



where R is 24 h respiration (kg CO2 m -2 d _1 ), P is daily gross 
photosynthesis (kg CO2 m -2 d _1 ), W is plant mass (kg CO2 
equivalents), k is a constant (dimensionless), and c is has 
units of d _1 . McCree applied a 12 h light period during a 
24 h day. Details of how R and P are obtained from the 
measurements are given by McCree (1970), a procedure 
which could be affected by daylength. McCree' s estimate of 



the parameters was k = 0-25 and c ■ 
through by dry mass W gives 



0 015 d \ Dividing 



R P 

— = k—- 
W W 



c= 0-25 — + 0-015 
W 



(2) 



This equation relates a specific 24 h respiration rate to a 
specific daily gross photosynthetic rate. The simplicity and 
intuitiveness of these equations has contributed greatly to the 
appeal of the GMRP paradigm, with c being interpreted as a 
maintenance coefficient. 



The situation here may be compared with micro-organisms 
in a chemostat (Pirt, 1965), where it is possible, in a true 
steady state without diurnality, to measure directly the gross 
rate of supply of substrate, respiration rate and also growth 
rate. This leads to a similar equation. 

The meaning of eqns (1) and (2) has been much discussed 
over the years (Amthor, 2000). A common interpretation is 
that respiration is viewed as having two components: one 
results from growth of the organism [growth respiration 
(GR)]; the other (MR) is attributed to the organism maintain- 
ing its status quo (Pirt, 1965; Thomley, 1970). 

Two objections can be raised against GMRP. The first is on 
scientific grounds: the biochemical and other processes occur- 
ring in a growing organism are qualitatively the same as those 
which occur in an organism which is maintaining the status 
quo. This especially applies to the metabolic processes 
which give rise to respiration and yield ATP and reducing 
power, but also those anabolic processes which are synthesiz- 
ing compounds which in some cases are replacing compounds 
that have recently been catabolized. In general, there is not a 
distinct set of processes which can be said to belong to 'main- 
tenance', and another set which belongs to 'growth', although 
catabolic processes (such as protein breakdown) may lead to a 
growth-associated process (protein synthesis) maintaining the 
status quo of the plant and giving a respiratory flux which 
can be dubbed maintenance. 

The second objection is a practical matter: many authors have 
used the paradigm when constructing models of plant growth, 
crop growth and plant ecosystems. Again and again unaccepta- 
ble behaviour of the model has occurred, and the maintenance 
component has to be arbitrarily 'fixed'. For example, de Wit 
et al. (1970, p. 61) remarked 'Simulation of this viewpoint 
leads to inconsistent results. If it is assumed that the respiration 
per unit plant material is low, it appears that yield levels which 
are observed in the field may be obtained, but then simulated res- 
piration rates under controlled conditions are far too small. If it is 
assumed that respiration per unit of plant material is higher, 
simulated respiration rates under controlled conditions may be 
in the observed range, but then simulated ceiling yields in the 
field are far too small.' Loomis (1970, p. 140, 141) states 'The 
McCree equation R = kP + cW summarizes well certain data 
. . . but it would be a mistake to assume it holds a priori for 
other situations, . . .', and 'The big difficulty is, of course, in 
accounting for 'maintenance' respiration...'. The impasse is 



Photosynthesis, P q 



Substrate 
mass 
M cs 



U, utilization = /c q M c 



Respiration, R= (1 - Y G )U 



Growth, 
►— 



Gy= Y n U 



Structural 
mass 



Senescence/ 
degradation, 
S=k x M cx 



Litter input, 
'ut=0-W)S 



Fig. 1 . Plant growth model for respiration. The two state variables are substrate mass, M cs , and structural mass, M cx - There are four key parameters: Y G , growth 
efficiency - fraction of substrate C utilized for growth which appears in structure - the remaining fraction (1 - Y a ) is respired; / rec , fraction of senescing/ 
degrading structural C recycled to substrate, the remainder becoming litter; k s , rate constant for substrate utilization; and k x senescence rate constant. 
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resolved (typically) by making the maintenance coefficient vari- 
able, e.g. Seginer (2003) assumes that the maintenance coeffi- 
cient depends asymptotically on plant structural mass. Van 
Oijen et al. (2010) assume that MR equals GR, an assumption 
tantamount to jettisoning GMRP in favour of a growth-only res- 
piration model. Others assume that the MR rate coefficient 
depends on the availability of substrates (e.g. Thornley, 1998, 
p. 38), as in the growth process. The problems revolve around 
issues such as: a common substrate (carbohydrate) is used for 
both growth and maintenance; priorities (more explicitly, utiliz- 
ation rate constants) have to be assigned between growth and 
maintenance; and these priorities need to change according to 
substrate supply and utilization (which determine substrate frac- 
tion). There is no straightforward way of dealing with growth 
and maintenance separately because the pools and anabolic pro- 
cesses are the same for both growth and maintenance. Many 
modellers do not like to represent substrates because their rep- 
resentation is perceived to be difficult (which it can be), 
although consideration of the differences between sub-arctic 
and semi-tropical plant ecosystems, and the diurnality of 
many plant variables including substrate pools might suggest 
that their representation is essential for realism. 

Here a two-state variable model (structure and substrate; 
Warren Wilson, 1967) is described. 'Substrate' includes 
storage (easily remobilizable) components. The model is a 
simplification and extension of a model proposed some years 
ago (Thornley, 1977; see also Loehle, 1982, who expanded 
on the 1977 model and provided a reconciliation of it with 
the more traditional view of respiration). The simplification 
is that non-degradable structure is omitted so that steady-state 
solutions can be derived, providing the insights that can be 
obtained from analytical expressions. The extension is that 
senescence is included, again so that steady-state solutions 
can be obtained. The significant processes are photosynthesis, 
growth with GR, senescence and recycling to the substrate 
pool. Respiration is an output of the growth process alone. 
Maintenance is not a feature. The model is simple enough 
(with three significant constants) to permit a thorough explora- 
tion of its properties, for both exponential growth and the 
steady state. My first objective is to demonstrate that its predic- 
tions are reasonably compatible with McCree's equation [our 
eqn (1)] (McCree, 1970), and therefore provides an alternative 
to the GMRP. My second objective is to define and extract a 
specific MR rate for a non-growing plant (which can occur 
in two ways). My last objective is to show that the model 
gives explicit expressions and acceptable and conservative 
predictions for the respiratiomphotosynthesis ratio (r R:/ > ), 
which for plant ecosystems is often in the range 0-4-0-5. 
Approximate expressions for the parameters of McCree's 
equation can be derived in terms of the significant parameters 
of the current model. 



MODEL AND METHODS 

The scheme is drawn in Fig. 1. Variables and parameters are 
listed in Table 1. The two state variables (pools) are Mqs 
(kg substrate C m -2 ) and M cx (kg structural C m -2 ); these 
are per unit ground area. Next the fluxes in Fig. 1 are calcu- 
lated, all having units of kg C m -2 (ground) d _1 . 



Photosynthesis, V g 

The expression used for gross photosynthesis [eqn (4)] has a 
long pedigree, and is parameterized with familiar quantities. It 
gives the correct dependencies on plant size (leaf area index), 
light flux density, photosynthetic efficiency and light-saturated 
leaf photosynthesis. These are needed because we shall be 
comparing our results with measured data. 

First calculate the leaf area index, Lai, from structural mass, 
M cx , with 



Lai = LarMqx- L 



AR 



50 m 2 leaf (kg structural C) -1 (3) 



Leaf area ratio, L AR , is a rounded value obtained from Causton 
and Venus (1981, p. 46, fig. 2-6a) who give a leaf area ratio of 
about 0-02 m 2 (g total mass) -1 . Note that L AR is a 'structural' 
leaf area ratio, with reference to just the structural biomass C 
component. Therefore, it is higher than usual, which is per 
unit of total biomass. We assumed that 40 % of total dry 
mass is structural C (50 = 1000 x 0-02/0-4). It is assumed 
that L AR is constant in order to make a simple analysis feas- 
ible. However, Bertin and Gary (1998) examined leaf mass 
per area (LMA, g m -2 ) in tomato, manipulating light, CO2 
and fruit load. They found substantial variations of incremental 
LMA related to source/sink status and sugar and starch 
content, and used their findings to modify their modelling of 
tomato plant growth. Such modelling requires keeping track 
of leaf area age distribution and increases complexity greatly 
(e.g. Thornley, 1998, pp. 48-49). 

Gross photosynthetic rate, P g (kg C m -2 d -1 ), is calculated 
with 



P g = ^86400[l-exp(-L AI )] 
12. 



akPn 



ak + /"max 



= 77 86400[1 - exp(-L AR M cx )] . 

44 al 0 + P max 

a= 10 x 10~ 8 kgCO 2 (JPAR)- 1 = 0 05molCO 2 

(quantum PAR) -1 

7 0 = 50 J (PAR) m" 2 s~ 1 = 230 yumol PAR irT 2 s" 1 

P max = 0-5x 10- 6 kgCO 2 irrV = ll-4/LunolC0 2 nr 2 s -1 

(4) 



P g is, via eqn (3), a function of the structure state variable, 
M C x- This equation is based on a continuous canopy of leaf 
area index L A i, light flux density incident on the top of the 
canopy of 7 0 , exponential decrease in light level on descending 
into the canopy of exp(- L AI ), a rectangular hyperbola for leaf 
response to light [P lea[ = aI 0 P ma J(aI 0 + P max ), kg C0 2 m" 2 
(leaf) s ] with photosynthetic efficiency of a, a light- 
saturated photosynthetic rate of Pmax at the top of the canopy 
which decreases as exp(-L A i) on descending in the canopy, 
and black leaves with a canopy extinction coefficient of 
unity [see, for example, Thornley and France, 2007, 
pp. 297-298, equation (9-50) with k=\ and m=0]. The 
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Table 1. State variables, parameters and definitions: relevant equation numbers are indicated 



State variable 


Description 


Initial (f = 0) value 


M C s 
Max 


Mass of substrate C [eqns (17), (1 1)] 
Mass of structural C [eqns (17)] 


0-000005 kg substrate C nT 2 
[= 0 05 M cx (t = 0)] 
0-0001 kg structural C nT 2 


(b) 


Parameter + value 


Description 


Units [with alternative where 
applicable] 


,/di = 0-5 


Fractional day length [(43)] 


h daylight (24 hr' 



/rec = 0-5 Fraction of senescence flux recycled [eqns (16), Fig. 1] 

lo = 50 [230] (continuous light) 100 Light flux density at top of canopy [eqns (4), (44)] 
[460] (12 h day) 

Rate constant for substrate utilization [eqns (12)] 
Senescence rate constant [eqns (15)] 
Leaf area ratio [eqns (3)] 

Saturating photosynthetic rate at top of canopy [eqns (4)] 



k x = 0-04 
L AR = 50 

P max =0-5 x 10 6 [11-4] 



F G =0-75 
a= 10~ 8 [0-05] 
/3 = 0-39 
p m = 00 127 



Growth efficiency [eqns (13)] 
Photosynthetic efficiency [eqns (4)] 

Constant for linearized gross photosynthetic rate at I Q = 100 [eqns (8)] 
Specific maintenance respiration rate [eqns (31), (42)] 



J (PAR) m 2 s 1 [amol (PAR) 

m" 2 s" 1 ] 

d" 1 

d-' 

m 2 leaf (kg structural C) 1 
kg C0 2 m" 2 s -1 [p.mol C0 2 
m" 2 s^ 1 ] 

kg structural C (kg substrate 

cr 1 

kg C0 2 (J PAR) -1 [mol C0 2 

(quantum PAR) -1 ] 

kg C substrate (kg structural 

cr 1 d- 1 



(c) 

Variable 



Description 



Units 



C s , C s (eg), C s (ss) 

/Rnight:R24h 

G x 

Ave 

Lai 
M c 
P,, P 



R, R„ 



P g (ss), P n (ss), 
„ fi(ss), R 24 h 



g,day? * n.day 



rR:P t , r R:Pi (eg), r R:Pi (ss), 

''R24h:Pgday 



' RnightiPnday. 



s 
u 

(J-s, n-x; 



|J-(eg) 

^g.day 



'n.day 



P. Pm. Pnighb PlAh 



Substrate fraction [eqn (11)], in exponential growth [eqns (25), (27)], in steady state 
[eqns (42)] 

Fraction of 24 h respiration occurring at night [eqns (47)] 

Structural growth rate (gross) [eqns (13)] 

Litter flux [eqns (16), Fig. 1] 

Recycled flux to substrate pool [eqns (16), Fig. 1] 

Leaf area index [eqns (3)] 

Total mass of system [eqns (18)] 

Gross, net photosynthetic rates [eqns (4), (19)]; steady-state values [eqns (36), (38), 
(41)]; daily values [eqns (45), (47)] 

Respiration rate [eqns (13)], night-time integral [eqns (47)], steady-state value 
[eqns (41)], daily value [eqns (45)] 

Ratios of: respiration to gross photosynthesis [eqns (21)], in exponential growth 
[eqn (34)], in steady state [eqns (42)], night respiration to daytime net 
photosynthesis [eqns (47)], 24 h respiration to daytime gross photosynthesis 
[eqns (46)] 

Senescence rate [eqns (15), Fig. 1] 
Utilization rate of substrate [eqns (12), Fig. 1] 

Proportional rates of growth [eqns (20)]; in exponential growth [eqns (25), (27)] 
Specific gross photosynthetic rate [eqn (6)], in exponential growth [eqn (33)]; daily 
gross, net values [eqns (48)] 

Specific respiration rate [eqn (14)], at maintenance [eqns (31) or (42)], night value 
[eqns (48)], 24 h value [eqns (48)] 



kg substrate C (kg structural 

cr 1 

kg structural C m~ 2 d~' 

kgCnT 2 d~' 

kg Cm" 2 d~' 

m 2 leaf (m 2 ground) 1 

kg total C irT 2 

kgCnT 2 d~' 



kgCnT 2 d~' 



kg structural Cm d 
kg substrate Cm^f 1 
d-' 

kg C (kg total Cr 1 d~' 

kg C respired (kg total C)~ 
d" 1 



factor 86 400 s d _1 converts per second units to per day; 12/44 
converts mass of CO2 to mass of C. The light flux density I 0 of 
50 here is for continuous light (giving a daily light receipt 
of 4-32 MJ photosynthetically active radiation (PAR) m -2 
d _1 or 19-9 mol PAR m -2 d _1 ); elsewhere a value of 100 is 
applied for a 12 h day, giving the same daily light receipt. 



Equation (4) is used because of its simplicity and tractabil- 
ity. There are other better approaches available (e.g. Thornley 
and France, 2007, pp. 288-301; Johnson et ai, 2010). 
Although plant growth does of course depend on the par- 
ameters in eqn (4), the key results of our analysis are indepen- 
dent of the parameters in eqn (4) and its precise form, as long 
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as the equation applied to calculate gross photosynthetic rate 
can be linearized for low values of L A1 so that exponential 
growth is a possible solution. 

It will be useful later when points of inflexion are calculated 
to have an expression of the differential of P g with respect to 
M cx , namely 



dP g 
AM, 



cx 



44 



-86400 



' ado + Pn 



-Lar exp(-L A RM cx ) (5) 



exponentially at a proportional growth rate dependent on I 0 . 
Later, we use the inverse of the second of eqns (8), namely 



/o 



P 



( 86400 ^W„ 



50 



P 



Anax " P 



50 



p 



0-589-/3 



(10) 



This decreases exponentially with M cx from 0-29 kg substrate 
C (kg structural C)" 1 d" 1 when M cx is small (7 0 = 50) at 
time t = 0 (Table la) towards zero as structural mass M cx 
(and L A1 ) — >• oo. 

The specific instantaneous gross photosynthetic rate, ir g [kg 
C (kg total C)" 1 d" 1 ], is 



1 



M c M cx + M cs M cx 1 + C s 



(6) 



Total mass Mc is given in eqns (18); see eqn (11) for substrate 
fraction Cs- Note that increasing Cs decreases the specific 
gross photosynthetic rate. 

At low and high values of light flux density I 0 , eqn (4) 
approximates to 



P g (/„ -» 0) = ^86400 o/ 0 [l - exp( - L AI )] 
P g (I Q -> oo) = ^86400 P max [l - exp( - L M )] 



(7) 



For low values of leaf area index, L M [or structural mass, M cx , 
eqn (3)], we write [from eqn (4)] 



P g = g 86400 aI ° P ™ 

g 44 Od 0 + ^max 

«/„) = 77 86400 ° mdx L AR 
44 al 0 + P max 



LarMqx = PMcx 



(8) 



to calculate light flux density I 0 
Decreasing j6 decreases I 0 . 



from a given value of /3. 



Substrate utilization, U, respiration, R, and structural growth, G x 

Substrate 'concentration' or fraction, Cs [(kg substrate C) 
(kg structural C) _1 ], is a useful intensive variable, defined as 



C s = M cs I Mcx 



(11) 



Assuming that the volume of the reaction in which substrate 
utilization takes place is proportional to the mass of structural 
C, M cx , then the rate of substrate C utilization, U (kg substrate 
C m -2 d -1 ), is (assuming linear or mass-action kinetics with 
rate constant, k s ) 



U = k s M C xCs = k s Mcs- k s = 1 d 



-l 



(12) 



With all defaults and in continuous light (Table lb), this k s 
value gives a substrate fraction Cs [eqn (11)] = 0-27 
during the initial exponential growth phase [eqn (25)], and 
C s = 0 053 at the steady state [eqn (42)]. A total of 63 % of 
an initial quantity of substrate is utilized within one 24 h day 
[exp(-* s ) = 0-37] [eqn (12)]. 

Utilized flux of substrate C, U, is divided between structural 
growth, G x , and respiration, R (kg C m -2 d _1 ) (Fig. 1) [note: 
we are not now using McCree's units for R; see eqn (1)]: 



G x =Y C U = Y G k s M cs and R = (1 - Y G )U 
= (1 - Y c )k s M C s- 7 G = 0-75 



(13) 



Parameter j8 has units of per day. In carbon units, it is the frac- 
tion of current structural C added per unit time via gross photo- 
synthesis to the C substrate pool. It is a rectangular hyperbolic 
function of light flux density I 0 . Values in the dark (7 0 = 0), at 
Io = 50 J PAR m -2 s _1 and its maximum are 



= 0) = 0 
/3(/ 0 = 50) = 0-29 



Pn 



m 



12 



OO) = 86400 — P ma x£AR 



0-589 d" 



(9) 



At high continuous light and low L AI , gross photosynthesis 
adds C equal to approx. 60 % of the current structural C to 
the substrate C pool per day. 

The linear expression P g = f}M C x of eqn (8) is valuable 
because it means that at low values of L A1 the plant grows 



Y G is a fraction (0 < Y a < 1) termed growth efficiency or 
growth yield. Y G values can vary considerably. The value of 
Yq = 0-75 is fairly typical for vegetative plant material (e.g. 
Thornley and Cannell, 2000; Thornley and Johnson, 2000, 
pp. 350-353). No respiratory flux is associated with the senes- 
cence processes in Fig. 1 . Such fluxes may exist and could be 
linked to outputs from a litter pool or to the inputs to that pool; 
they are omitted so that we retain the simplicity of a single res- 
piratory output flux. 

The instantaneous specific respiration rate, p [kg C respired 
(kg total C mass)" 1 d" 1 ], is [cf. eqn (6)] 



(1 - Fg)^sMcs (1 - Y G )k s C s 



R 



M c M cs + Mcx 



1 + Cs 



(14) 



For total mass M c see eqns (18). p increases from zero to a 
maximum of (1 - Y G )k s as Cs increases from zero to oo. 
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Senescence flux S, fluxes to litter, l Ut , and recycling of substrate, l rec (kg C m 2 d : ); therefore, 

The proportional rate at which structure senesces is k x 
(d _1 ), giving a senescence flux of S (kg structural C m -2 



d"\ Fig. 1) of 



S = k x M CX - k x = om<r 



(15) 



The fraction senescing or degrading per day is 4 %, giving a 
tissue lifetime of 25 d. No distinction is made between senes- 
cence and (non-pathological) degradation, and henceforth the 
process is referred to as senescence. More detailed plant 
models may distinguish between ageing and senescence, and 
non-age-related tissue turnover or degradation (e.g. Thornley, 
1998). No respiratory flux is associated with senescence (see 
above). The structural C flux from senescence is divided 
between a fraction, / rec , which is recycled to the substrate C 
pool with flux, I rec , and the remainder, 7 lit , which is input to 
a litter pool (Fig. 1). Thus 



M c = M cs + M cx 
dM c dMcs dM C x 



At 



dM c 
~dT 
dM c 
~dT 



Pg(M C x) ~ k s M cs +f rec k x M cx 



dt ' dt ~ g 
+ Y c k s M cs - k x M C x 

PgWcx) -k s (l- Y G )M CS - k x (l -f iec )M cx 

P g (Mcx) -R- /lit = Pn - /lit, where^ = P g -R 

(18) 



Equations (13) and (16) have been applied. Net photosynthetic 
rate, P n , is defined as gross photosynthesis P g [eqn (4)] less 
respiration R [eqns (13)]. These two alternative forms for net 
photosynthesis, namely 



P n =P g -R = 



dMr 



dt 



(19) 



= frecS =f rec k x M cx and I lit = (1 -/ rec )S 
= (1 -f KC )k x M cx . = 0-5 



(16) 



It is assumed that a constant 50 % of the senescing C is 
recycled. In general, it may be expected that this is a variable 
fraction depending on nutrient status (level of substrate frac- 
tion Cs) and tissue composition. Note that 0 < f rec < 1 . In a 
real varying environment where substrate supply moves 
between scarce and plentiful, it does not seem unreasonable 
to take = 0-5 as representing an average value. 



Differential equations 

By inspection of Fig. 1 and using eqns (4), (12), (13), (15) 
and (16), these are 



dM cs 
dt 



Pg(M C x) - k s M cs +f KC k x M C x, 



M cs (f = 0) = 0-000005 kg substrate C m" 
dM cx 



dt 



= Y c k s M C s - k x M cx , 



M cs (t = 0) = 0-0001 kg structural C 



(17) 



Gross photosynthetic rate P g is written as P g (M cx ) to empha- 
size that P g is a function of state variable M cx [Fig. 1, eqns (4), 
(3)]. The t = 0 value of substrate fraction Cs is 5 % [eqn (11), 
defining M cs (?=0) in terms of M cx (f = 0)]. Denote total 
mass by Mc (kg C m~ 2 ) and its rate of change by dMc/df 



are used to check model formulation (stoichiometric validity) 
and programming. 

Proportional growth rates of state variables (jul s , |jl x ) and 
total mass (|x)(d~ ) are 



Ms 



1 dM cs 
M cs dt 

1 dM cx 
' M cx dt 
1 dM c 



M c dt 



(20) 



These are intensive variables which are constant during 
exponential growth. 

The instantaneous ratio (r) of respiration R [eqns (13)] to 
gross photosynthesis P g [eqn (4)], r R:P , is 



R 

r R: p g - y 



(21) 



The equation for the double differential of total system C is 
[differentiating the third of eqns (18)] 



d 2 M c dP g dM cx 



dt 2 



dMcx df 
-(1 -/«)** 



(1 - Y G )k s 

dMcx 
dt 



dM c& 
dt 



(22) 



With eqns (17) and (5), d 2 Mc/df 2 can be written in terms of 
the state variables. Hence time of inflexion when d 2 M c / 
df 2 = 0 and state variable values at inflexion can be extracted 
numerically. 

The differential eqns (17) could be non-dimensionalized 
[using eqn (4)], but not uniquely. This was not seen as a 
helpful simplification. 
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Exponential growth (eg) 

There are two limiting cases for which analytical solutions 
exist. These are first exponential growth and secondly the 
steady state in a constant environment. In this section we 
examine exponential growth (eg). Exponential growth occurs 
at low L A1 [eqn (8)] when eqns (17) become 



dM cs 

dt 
dMcx 

dt 



= P(Io)Mcx - k s M c & +f KC k x M cx 
= Y G k s M c& - k x M cx 



(23) 



/3(/o) is so written to make explicit that p is a function of light 
flux density I 0 [eqns (8)]. Assuming constant and equal 
proportional growth rates |x for M C s and M C x [eqns (20)] 
and constant substrate C fraction, Cs [eqn (11)], eqns (23) 
become (divide through by M cx ) 



lxC s : 
H 



P(k) + &x/rec 

Y a k s C$ - k x 



ksC s 



(24) 



These equations are solved for |x and C s . Eliminating jjl (the 
positive square root gives the biologically meaningful result) 



Y G k S C| + (k S - k X )C S - [j8(/ 0 ) + Wrec] = 0 

a = Y G k s , b = (k s - k x ), c = -[(3(I 0 ) + kxf rec ] 



Cs(eg) 



2a \ 



4ac\ 



/i(eg) = Y G k s C 5 (eg) - 



(25) 



Here, 'eg' denotes exponential growth. 

Plotting exponential growth rate |Ji(eg) against substrate 
fraction C s (eg) gives a straight line [last of eqns (25)], with 
zero growth occurring when 

F G A: s C s (eg; fi = 0) = k x ; C s (eg; n = 0) = k x /Y G k s 

= 0-05333* (26) 

There is a minimum (non-zero) substrate fraction of 0-05333* 
kg substrate C (kg structural C) _1 (default parameters, 
Table lb) below which the plant dies [|x(eg) < 0]. Low 
values of the senescence rate constant k x and high values of 
growth efficiency Y a and substrate utilization rate constant k s 
lead to a low value of Cs(eg; |x = 0) [see eqn (30) below]. 

In terms of the underlying parameters, Cs(eg) and |Ji(eg) 
become [eqns (25)] 



Cs(eg) 



1 



2Y G k s 



- k x ) 2 + WMWo) + kxf iec ] - (As - k x )) 



Meg) 



(27) 



(VWs - k x f + 4Y G k s [(3(I Q ) + kxfrec] - (h + k x )) 



For small values of (3(Io) [which is proportional to light flux 
density, I 0 , eqn (8)], these are linearly proportional to f3 and 



thence 7 0 [because ^/(l + x) ~ 1 + for small x]. For large 
/3, Cs(eg) and |x(eg) ~ ~Jf3. Remember that /3 has an asympto- 
tic hyperbolic dependence on light flux density 7 0 , with asymp- 
tote /3 max depending on maximum photosynthetic rate P max 
[eqns (8), (9)]. 

In the dark (7 0 = 0), /3 = 0 [eqns (8)], eqns (27) become 



C s (eg;/ 0 = 0) 



21gA:s 

(V(*S - kxf + 4Y G k S kxfrec - (*s - k x )) 

= 0-0205 kg C substrate (kg structural C)" 1 



(28) 



M(eg;/o= 0) = - ( V(*S - kxY + 4y G *s*x/n* - (*s + kx) 
= -0-0246d _1 



Assuming that k s ^> k x (Table lb), these can be approximated 
by 

C s (eg;/ 0 = 0;£ s » k x ) » f lec k x = 0 -02 kg C substrate 

(kg structural C) _ 1 
Meg;/,, = 0;k s » k x ) « k x (l -f KC Y G ) = -0-025 d" 1 

(29) 



The rate of decline in the dark is least if growth yield Y G and 
recycling fraction f Iec are high (approaching unity). 

The plant dies if jjl < 0. For (3 the second of eqns (27) gives 
an upper limit, and then eqns (10) and (9) give an upper limit 
of light flux density 7 0 , so that |x < 0 implies 

(k S - k x ) 2 + 4Y G k S [P + fcx/rec] < (*S + kxf 

P < — ^2 = 0-03333 



/o <^ 50 P 



50- 



0-03333 



a P max -P p mwL (= 0-589) - 0-03333 
= 3 J PARm" 2 s _1 [14 x 10 _6 molPARm" 2 s _1 ] 



(30) 



[This result can be derived directly from eqns (24) (eliminate 
Cs; write down the condition for (jl < 0).] Death of the plant 
depends on light flux density I Q [eqn (8)], but the critical 
values of /3 and I 0 for positive growth are less (giving shade 
tolerance) for lower values of senescence rate parameter k x , 
and for higher values of growth yield Y G and recycling fraction 
f mc . There is no dependence on substrate utilization rate 
constant k s [eqn (12)]. 

At zero growth rate, the plant is just maintaining its status 
quo. The specific maintenance respiration rate [eqn (14)], 
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p m , is [with eqn (26)] 



p(eg;p= 0) = p, 
= 0-0127 cT 1 

or for ^ s » k x , p(eg; p. = 0) 
= 0-0133* cT 1 



R k s k x 
— = (1 - Y G ) — — 

m c u; Y G k s + k x 



l-Y G 
Yc 



(31) 



Values for default parameters (Table lb) are given. There is 
weak dependence on utilization rate k s [eqns (12)] as generally 
Yaks ^> kx- Pm does not depend on recycling fraction / rec 
(Fig. 1) [but see Fig. 4B and eqns (56) - the intercept of 
McCree's equation does depend on the recycling fraction / rec ]. 

The specific respiration rate for a plant in exponential 
decline in the dark is [eqns (14), (28)] 



p(eg;/ 0 = 0) = 



R (l-Y G )k s C s (eg;I 0 = 0) 



M, 



c 



l + C s (eg;/ 0 = 0) 



0-00502 d" 



(32) 



For k s ^> k x , this is approximately (1 - Y a )f Tec k x = 0-005 d~ 
[use eqns (29)]. This quantity is sometimes assumed to be a 
measure of the MR rate, but it is substantially less than the 
specific respiration rate required to maintain the plant's 
status quo [eqns (31)]. 

The specific gross photosynthetic rate, 7r g [eqns (6), (8), 
(13)], for exponential growth, is 



^g(eg) 



/3(/o) 



1 + C s (eg) 



(33) 



For a given light level and value of /3 [eqn (8)], ir g decreases as 
the sugar level in the plant [Cs(eg), eqn (27)] increases (due 
perhaps to a decreasing A: s ; Fig. 1, eqns (12)]. 

The ratio of respiration to gross photosynthesis (19) is 
[eqns (13), (8), then (27)] 

R (l-Y G )k s M cs (1 - Y G )k s C s (cg) 

>R:P E (eg) = — (eg) = 



PMcx 



= (1 - Y C ) 
2Y G P 

= 0-23 



(k s - k x f + 4Y G k s (l3 + kxf rec ) - (k s - k x ) 



(34) 



given is for Iq = 50 W 
s" 1 ) giving (5 [eqn (8)] = 
At zero growth rate this becomes 



The numerical value 
(230 iJimol PAR m" 2 



(PAR) m" 
0-2945 d" 



ni:P E (eg;p = 0) 



1-Yq 
1 - Ycf m 



= 0-4 



(35) 



Here we used the second of eqns (30) with the equality to sub- 
stitute for j3 and eqn (26) for Cs in the first line of eqns (34). 
The ratio at zero growth rate (maintenance) is 0-4 with default 
parameters (Table lb) or 0-25 with no recycling (/ rec = 0) [see 
paragraph after eqns (16); also eqns (30) and following 
paragraph]. 



Steady state (ss) 

In eqns (17), put the derivatives equal to zero, giving 
(ss denotes the steady state) 



0 = Pg(ss) - £ s M cs (ss) +/ rec &xMcx(ss) 
0 = Y c k s M C s(&&) - ^xMcx(ss) 



(36) 



Solving for M cs , M C x an d adding to give total mass M c , 
therefore 



Mcs(ss) 
\'h (ss) = 



P g (ss) 



■ , Mc X (ss) 



*S(1 -frecY G ) 

Mcs(ss) + Mcx(ss) = 



y G P g (ss) 

feed -/rec^G) 

(k x + k s Y G )P ^ g (ss) 
^s^x(l -f K cY G ) 



(37) 



To find the steady state, solve two equations for -P g (ss) and 
M cx (ss) [from eqn (4) and the second of (37)]: 



P g (ss) = ^86400(1 - exp{ - [L AR M cx (ss)]}) 



44 

n , , ^x(l — frecY G ) 

P s (ss) = Mcx(ss) 



ado 



(38) 



Eliminating P g (ss) gives 

k X (l -/rec^G) 



12. 



^G 



Mc X (ss) = —86400 
44 



(39) 



(1 - exp{ - [L AR Mcx(ss)]}) 



al 0 + P n 



This transcendental equation can be solved numerically for 
M cx (ss), e.g. using the Newton-Raphson method (Thornley 
and France, 2007, pp. 830-831), noting that eqn (39) is 
easily differentiated with respect to M C x(ss). There is a 
trivial solution at M C x = 0. There is a non-zero solution 
only if, at small M C x, the slope of the right side of eqn (39) 
is greater than that of the left side. Therefore, using eqn (8) 
(P g = jSMcx), 



&X(1 -/recic) 

/3 > — = 0-03333 



(40) 



This inequality is equivalent to eqn (30) for plant death during 
the exponential growth phase. With eqn (10), it is equivalent to 
/ 0 (jjl > 0) >3 J (PAR) m" 2 s" 1 [14 x 10" 6 mol (PAR) m" 2 
s _1 ] for standard parameters (Table lb) [cf. eqn (30) and 
following paragraph]. 

Given that eqn (40) is satisfied and a non-trivial steady-state 
solution to eqn (39) exists for M cx (ss) and therefore also 
P g (ss) with either of eqns (38), the respiration rate R in the 
steady state becomes [eqns (13) and (37)] 



R(ss) = £ s M cs (ss)(l - Y G ) 



l-Yc 

1 —fiecY G 



P g (ss) (41) 
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Hence follow several important ratios for steady- state plant 
ecosystems [for Cs(ss), use eqn (11) and the second of 
eqns (36)], with parameters from Table lb: 



C s (ss) 



Mcs(m) 



M cx (ss) 

rR:P g (SS) = — (SS) 



p(ss) 



R(ss) 



kx 

^- = 0-05333 
Y G ks 

- 1 ~ YC = 0-4 

1 —frecYa 

k s k x (l - Y a ) 



M c (ss) k x + k s Y G 



0-0127 d" 



(42) 



These are all independent of environment, here light flux 
density, in contrast to these quantities during exponential 
growth [eqns (27), (34)] where they depend on /3 and light 
flux density I 0 [eqn (8)] apart from the specific MR respiration 
rate, p m [eqn (31)]. The last of eqns (42), for specific respir- 
ation rate in the steady state, is obtained from the first of 
eqns (41) and the last of eqns (37). This is identical to eqn 
(31) for the specific respiration rate of a plant growing expo- 
nentially at zero rate. This important general result provides 
a link between long-term ecosystem behaviour and short-term 
experimentation on young plants. The R'-P g ratio in eqns (42) 
has a simple interpretation with the Fig. 1 model because the 
denominator, 1 - x (say), can be expanded as 1 + x + x 2 + 
x 3 . . ., giving successive cycles of substrate use and re-use. 



gross photosynthesis is 



^24h 

;,day 



(46) 



These quantities are calculated so that a comparison (albeit 
imperfect) with McCree's (1970) results can be made. 

CO2 gas exchange measurements generally give directly net 
photosynthesis during the day and respiration during the night. 
Again for comparison, and using an obvious notation, eqns 
(45) and (46) can be extended (including the fraction of 24 h 
respiration occurring during the night) [eqns (19), (13)] 



0.5 



n.day 



0 

_ ^night 

7"Rnight:Pn,day — p , /Rnig 

*n,day 



J P B dt, /? night = J Rdt, 
0 0.5 



R 



ht:R24h = 



night 



R 



24h 



(47) 



Specific values of daily photosynthesis and respiration are cal- 
culated by dividing by total mass M G [eqns (18)] at the end of 
each 24 h day: 



77"n,day 
Pnight 



n.day 



^g.day = 



;,day 



M c 

^night 

~M~c' P24h ~Mc 



M c 

Rl4h 



(48) 



Diurnal environment 

Assume that light flux density 7 0 [eqn (4)] is varied on a 
daily basis, with 12 h of light followed by 12 h of darkness. 
The analysis of the last two sections dealing with exponential 
growth and a steady state in a constant environment is no 
longer possible. However, this more realistic situation can be 
examined numerically, obtaining results that may be compared 
with measurements such as those by McCree (1970) as well as 
with our analytical results for continuous light. 

Define fractional daylength, fa, by 



fa = 12 h daylight (24 h)" 1 = 0-5 



(43) 



Light flux density I 0 applied in eqn (4) to calculate photosyn- 
thesis is adjusted by dividing by fractional daylength: 



^0 := lo/fa 



(44) 



The symbol ': = ' means the value is replaced by. This ensures 
that daily light receipt is the same and makes simulations of a 
constant light regime and a 12 h day more comparable. 

A critical issue is how, for a diurnal environment, are 'daily 
gross photosynthesis' and 'daily respiration' defined. These 
quantities, denoted by P gi da y (kg C m -2 d _1 ) and respiration, 
^24h (kg C m" 2 d" 1 ), are calculated with [eqns (4), (13)] 



An exponential growth phase still occurs, now on a daily 
(24 h) basis rather than an instantaneous basis; there are fluc- 
tuations in proportional growth rates during the day 
[eqns (20)]. The steady state becomes a 24 h equilibrium 
state with diurnal variation. 



Re-parameterization of the model 

A more traditional single-state variable plant growth 
equation [cf. eqns (17) with two state variables] is obtained 
by making a zero-pool approximation (Thornley and France, 
2007) and by eliminating one parameter in favour of a syn- 
thetic (emergent) McCree-type parameter. 

Specific MR rate, p m (d _1 ), is [eqns (31) or (42)] 

k s k x (l - Y c ) 
Pm = p(eg; p = 0) = p(ss) = S X , , t (49) 

^G^S + Kx 

Any of the three parameters: k s , k x or Y G (Fig. 1) can be elimi- 
nated from the model in favour of p m . We eliminate k x , writing 



P m Yck s 



(1 - Y G )k s - p„ 
The differential eqns (17) become 



(50) 



g.day 



0-5 1 1 

f f f dMcs 

Pgdr = P g dr, R 24 h = Rdt (45) = -Pg(^cx) - k s M cs +f t 



PmYcks 



Daylength is 12 h and light is applied from 0 to 12 h and dark- 
ness from 12 to 24 h. The ratio of 24 h respiration to daytime 



dt 
dM cx 
dt 



(1 - Y G )k s - p„ 



-M, 



cx 



Y G k s M C s - 



PmYcks 



(1 - r G )* s - A 



-M cx 



(51) 
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For large values of k s , substrate content M cs is small, and is 
treated as a 'zero pool' so that inputs = outputs, dM cs /dr ~ 0. 
Also assuming that (1 - Y G )k s » p m [eqns (13), (12), (31), 
(42)], the first equation becomes 



(11-25 min). A program listing can be obtained from the 
author. 



RESULTS AND DISCUSSION 



0 = Pg(Mcx) - k S M CS +/rec 



Pm Y C 
(1 - Yq) 



Mcx (52) Standard growth curve 



Substituting for k$Mcs from eqn (52) in the second of eqns (51), 
and again for large k s [eqn (49) is used for the specific MR rate, 
Pm]: 

dMcx = Fofp g (Mcx)+/ rec ^^M cx V^M cx 



dt 

dMcx 
dt 

dMcx 
dt 

Finally 



(1-1g) 



i-y q 



Y G P g (Mcx) " (1 " Ycf zec )p^M C x 
1 — I G 



Y G ( P g (M cx ) - P m \ Y i" C M cx 



dM, 



cx 



dt 

m = p m Y G 



l-Y Q 



Y c P g (M C x) - mMcx 
1 - Yafn 



(53) 



l-Y G 



(1 - Ycf rec )kx 



(54) 



With default parameters and a constant environment 
(Table lb), plant mass grows as in Fig. 2A. This 
'expo-asymptotic' curve has a period of almost constant expo- 
nential growth from time f = 3 d to lOd when the L A1 reaches 
002. After a maximum, growth rate rapidly declines. 
Asymptotic L A1 (not shown, but see Fig. 6A with I 0 = 50) is 
8-8 (it is approx. 6 for a 12 h day). This was felt to be accep- 
table for a simple model simulating vegetative growth. In 
Fig. 2, the asymptote is at final mass Mc = 0T861 which is 
3-83 times the mass at the inflexion point. Logistic and 
Gompertz growth equations gives 2 and 2-7818 = e for this 
ratio (Thomley and France, 2007, p. 144, equation 5-19 and 
p. 147, equation 5-28). 

In Fig. 2B the ratio of respiration R [eqns (13)] to gross 
photosynthesis P g [eqns (4), (21)] is drawn. After an unimpor- 
tant transient, it increases from a value during exponential 
growth of 0-23 [eqns (34)] to an asymptotic value of 0-4 
[eqns (42)]. 



This is a single-state variable plant growth equation. The coeffi- 
cient (for k s = 100 d _1 ), m, is 



Proportional growth rate, (1/M C ) d/W c /df 



m = 0-025 d" 



(55) 



This can be compared with the specific MR rate of p m = 
001 33* d _1 [eqn (31)]. It is mistaken to interprets as the main- 
tenance requirement. Integration of eqn (54) or of eqns (17) with 
k s = 100 d _1 (say) gives exactly the same time course forM cx , 
so long as respiration is still defined by eqns (13). 

This result depends on the assumptions in Fig. 1. The 
crucial assumption concerns respiration rate R [here given by 
eqn (13) with k s M C s now given by eqn (52) in the zero-pool 
approximation]; the rest is stoichiometry. Equation (13) is 
what Warren Wilson et al. (1986) call a 'conditional' equation 
which represents a scientific hypothesis. When C conservation 
(stoichiometry) is considered, many equations that are 'identi- 
ties' can be written down which are algebraically self-evident 
[e.g. eqns (19)]. The terms 'conditional equations' and 'iden- 
tities' are sometimes referred to as 'synthetic statements' and 
'tautologies' (e.g. Ayer, 2001). Van Oijen et al.'s analysis 
(2010) consists mostly of identities, although two conditional 
equations are introduced: for respiration R (R = MR + GR), 
and then MR = GR. 



Numerical methods 

Programming was in ACSL (Advanced Continuous 
Simulation Language, Aegis Research, Huntsville, AL, USA; 
version 11-8-4), a fortran-based ordinary differential equation 
(ODE) solver. The ODEs are integrated using Euler's 
method with a fixed integration interval of At = 1/128 d 
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Fig. 2. Growth curve in a constant environment with default parameters 
(Table lb). (A) Plant mass, M c , eqn (18); growth rate, dM c /dt, eqn (18); and pro- 
portional growth rate of total mass, (1/M C ) dM c /dt, eqns (20). The black circle 
shows the inflexion point on the graph for M c (continuous line). (B) Ratio of res- 
piration R [eqns (13)] to gross photosynthesis P g [eqns (4)], r R:P [eqn (21)]. The 
dashed lines show the value of the ratio during exponential growth (0-23) and at 
steady state (04) [eqns (34), (42)]. 
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Exponential growth 

The analysis for an exponentially growing plant in a con- 
stant light environment is given in eqns (23)-(35). This corre- 
sponds to the 3-10 d time period shown in Fig. 2. In Fig. 3, 
light flux density I 0 , and thereby parameter /3 [eqn (8)], is 
varied across a range of I 0 from near zero to 500 J (PAR) 
m -2 s _1 (2300 p-mol m -2 s _1 ). Results are plotted against 
specific gross photosynthetic rate, ir g [eqns (6), (33)]. Two 
values for the substrate utilization constant A: s [Fig. 1; 
eqn (12)] are applied to demonstrate the effects of low and 
high substrate content (Fig. 3D). 

In Fig. 3A, specific growth rate |x [eqns (20)] increases with 
specific gross photosynthetic rate ir g , but at an increasing rate 
(Fig. 3E). This is because the substrate fraction [Cs, eqns (11), 
(27)] increases with ir % (Fig. 3D) and there is no yield factor 
such as Y G (Fig. 1) operating for the substrate component of 
dry matter. At low light (low 7r g ) and for high lc s , the substrate 



fraction is low; nearly all growth is structural which incurs the 
Y G penalty. The k s = 5 d _1 line in Fig. 3E approaches Yq = 
0-75 at low light (low 7r g ). 

Respiration is a decreasing fraction of increasing photosyn- 
thesis (Fig. 3B, C, F). The decrease in this fraction is marked 
for the k s = 1 d _1 lines because the substrate fraction 
increases to high values (Fig. 3D). The specific respiration 
rate p [Fig. 3B, eqn (14)] when growth rate |x is zero 
[eqn (31)] has values of 0-0127 and 0-0132 d" 1 for k s = 1 
and 5 d" 1 . 

Corresponding to the decreasing respiration fraction of 
Fig. 3B and F, the ratio of respiration to gross photosynthesis 
[eqn (21)] (Fig. 3C) decreases with increasing light. At zero 
growth rate this has a value of 0-4 [eqn (35)], decreasing for 
higher growth rates as the substrate fraction increases 
(Fig. 3D). The higher value of the substrate utilization constant 
k s = 5 d _1 gives lower values of the substrate fraction 
(Fig. 3D). Plant substrate fractions can cover a wide range, 
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Fig . 3 . Exponential growth solutions [eqns (23) to (35)] in constant conditions. Light flux density, 7 0 , and hence parameter /3 [eqn (8)], are varied from near zero 
to 500 J (PAR) m 2 s 1 (2300 (xmol PAR m 2 s '). Two values of parameter k s are applied. Specific gross photosynthetic rate ir g is given by eqn (33). The 
vertical dashed line labelled |ul = 0 marks the n g value for zero specific growth rate, jjl = 0, for k s = 1 d _1 (the value is slightly different for k s = 5 d ) 
[eqns (33), (30), (27)]. (A) Specific growth rate p., eqns (27); (B) specific respiration rate p, eqns (14), (27); (C) respiratiomgross photosynthesis ratio r R:P , 
eqns (21), (13), (8), (27); (D) substrate fraction Cs, eqn (27), (8); (E, F) slopes of graphs in (A) and (B), respectively, were calculated numerically. 
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Fig. 4. Simulation of McCree's equation [eqn (2); McCree, 1970]. A 12-h day- 
length is applied followed by 12-h darkness. In (A), (B) and (C), light during the 
day was varied from a low value to 1000 J (PAR) m~ 2 s _1 (4600 ixmolPARm -2 
s ). Two values of utilization rate constant k s were applied, giving high and low 
substrate plants. Plants are in 24 h exponential growth, attained after 10 d from 
low initial values (Table la). (A) McCree's formulation [eqns (1), (2); 
McCree, 1970] is followed. (B) Slopes of the lines in (A) are calculated numeri- 
cally. (C) The fraction of 24 h respiration occurring during the night is calculated 
assuming eqn (13) for respiration R is valid during the day and night 
[eqns (47)]. (D) Time courses of the substrate fraction are illustrated for a 
daytime light flux density of I 0 = 100 J (PAR) m~ 2 d -1 (460 ixmol PAR m~ 2 
d ). Specific 24-h respiration p24h an d specific daily gross photosynthesis 
7r g day are calculated in eqns (48). The vertical dashed lines labelled la = 0 in 
(A), (B) and (C) indicate where daily specific growth rate jjl = 0. 



from a few per cent up to 20-30 %, depending on species, 
growing conditions and time of day and season. 

McCree's experiment 

Now we apply the model to McCree's (1970) equation [our 
eqn (1)], attempting to follow his experimental protocol 
closely. No parameter tuning was attempted and default par- 
ameters are used (Table lb). In eqn (4), light flux density I 0 
is applied at a constant level for 12 h between the hours of 0 
and 12 h; otherwise between 12 h and 24 h, I Q = 0. It is 
assumed that McCree's plants were approximately in exponen- 
tial growth, and here the simulated plants are in 24 h exponen- 
tial growth. This is obtained accurately enough at time t = 10 d 
by using the initial values (/ = 0) of Table 1 a, and lasts until 
the leaf L Al is approx. 0 02, which depends on light flux 
density 7 0 . As in Fig. 3, two values of the substrate utilization 
rate constant A: s are applied (Fig. 1), so that the effects of high 
(with A: s = 1 d _1 ) and low (with A: s = 5 d _1 ) substrate fractions 
can be seen (Fig. 3D). Daily gross photosynthesis P a d ay (kg C 
m -2 d _1 ) and daily (24 h) respiration 7?24h (kg C d _1 ) are 
obtained with eqns (45). P g 4 ay and R 2 4h are both divided by 
the plant mass, M c [kg C m -2 , eqn (18)], at 24 h to give 
specific daily gross photosynthesis, tr g d ay (d _1 ) and specific 
daily respiration, p 2 4h (d _1 ) [eqns (48)]. It is assumed that res- 
piration R continues according to eqn (13) during both day and 
night (but see below). A range of daytime light flux densities 
7 0 [eqns (8)] from near zero up to 1000 J (PAR) m~ 2 s -1 
(4600 (xmol PAR m -2 s _1 ) gives a range of 7r g day . Low 
values of I 0 can give negative daily growth rates. 

Figure 4A shows respiration vs. photosynthesis for two values 
of A: s . The lines are curvilinear, but assuming linearity both have 
an intercept of approx. 0 005 d _1 . When slopes are plotted as in 
Fig. 4B, the lack of linearity is more apparent, especially for the 
lower lc s value which gives higher substrate fractions (Fig. 4D). 
The curvilinearity is less at high k s values where substrate frac- 
tions are smaller (Fig. 4D). Comparing Fig. 4A with McCree's 
equation in the form of eqn (2), our intercept of 0 005 d~ of a 
curve may be compared with McCree's value of 0 015 d _1 
from fitting a straight line to data (McCree, 1970, p. 227). 
Drawing a tangent to the k s = 1 curve in Fig. 4A at higher 
values of abscissa will give substantially higher values of the 
intercept. For low substrate (£ s = 5 d~ ) the slope decreases 
from approx. 0-248 to 0-245 (Fig. 4B) as photosynthesis 
increases; for high substrate (&s = 1 d _1 ) the slope decreases 
from 0-25 to 019. This occurs because growth in substrates 
does not incur a respiration penalty as does growth in structure. 
These slopes can be compared with McCree's estimate of 0-25 
(McCree, 1970, p. 227). As no parameter tuning has been 
applied, the general agreement with McCree's values is 
reasonable. 

Assuming that respiration throughout the 24 h is given by 
eqns (45) and (13) (but see below), a fraction for night-time 
respiration [eqns (47)] can be calculated (/? n i g ht/7?24h)- This is 
plotted in Fig. 4C. For the high substrate plants (& s = 1 d _1 ) 
this is roughly constant at about 0-5 (the small increase 
above 0-5 shown in Fig. 4C is due to the way in which the 
simulations were made: night follows day and the plants are 
larger (in structural terms) during the night due to growth 
and therefore respire slightly faster, although average substrate 
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levels at night are always lower than during the day [eqns (13) 
for R, Fig. 4D]. For low substrate plants (k s = 5 d _1 ), night 
respiration is well below the 0-5 expected from daylength 
alone, because for high k s average substrate level at night is 
much less than during the day (Fig. 4D). 



Sensitivity of McCree equation parameters to key model 
parameters 

Now we use the method of the last section (examining 24 h 
exponential growth solutions) to investigate the sensitivity of 
McCree equation parameters [McCree, 1970, p. 227; our eqn 
(2)] to key model parameters. These are: growth yield, Y G , 
recycled fraction, / rec ; senescence rate constant, k x \ utilization 
rate constant, k s ; and fractional daylength, / dl (Fig. 1). 

First from Fig. 4A and B, we note that McCree' s linear 
equation [eqn (2)] is only approximately satisfied by the 
model's predictions, which are in general curvilinear. 
Therefore, exact equations relating McCree parameters k and 
c of eqn (2) to model parameters (Fig. 1; Table lb) cannot 
be expected. 

Growth yield Y G (Fig. 5A) affects both slope and intercept 
[eqn (2), k, c] of specific 24 h respiration vs. specific daily 
photosynthesis [eqn (48)]. The recycling fraction f rec 
(Fig. 5B) affects the intercept proportionately but not the 
slope. With /j. ec = 0, the intercept is indeed zero. Senescence 
rate constant k x (Fig. 5C) similarly affects the intercept (pro- 
portionately) but not the slope. The utilization rate constant k s 
(Fig. 5D) has asymptotic effects on slope (see also Fig. 4A, B) 
and intercept (as A: s —> oo the lines superimpose). A k s value of 
50 d _1 gives a response that is little different from the k s = 1 
d _1 response. However, lower values of A: s (<1 d _1 ) cause 
substrate fraction Cs to be high (Figs. 3D, 4D); no respiration 
is associated with growth of substrate so this decreases specific 
respiration at higher values of specific gross photosynthesis; 
responses become more convex and the slope and intercept 
are decreased. In Fig. 5E, fractional daylength (/ dl ) is 
altered. Light flux density I 0 [eqn (4)] is adjusted so that the 
total daily light receipt is the same. i.e. light flux density I 0 
is twice as large for fa = 0-25 (a 6 h day) as for f dt = 0-5 (a 
12 h day). Because of the non-linear nature of eqn (4), daily 
gross photosynthesis is less for / dl = 0-25 than for f di = 0-5. 
Perhaps surprisingly, neither intercept nor slope is affected 
as daylength is altered. However if night-time respiration is 
plotted against daytime net photosynthesis [eqns (47)] 
then both slope and intercept are proportional to daylength 
[eqns (58)]. 

The results in Fig. 5 suggest that the two parameters of 
eqn (2) are given roughly by 



k « (1 - F G ), c « k x (l - Y G )f ie 



(56) 



Slope k depends only on yield constant Y G (Fig. 5A). Y G 
affects intercept c, which is also affected by senescence rate 
k x and recycling fraction / rec (Fig. 5B, C). In eqns (56) no 
dependence on the substrate utilization rate constant k s is 
included, which is reasonable if ks = 1 d _1 or higher. 

Equations (56) suggest that McCree' s k parameter [eqns (2)] 
can be interpreted directly in terms of growth efficiency Y G . 
Perhaps more interesting is McCree' s c parameter, often 
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Fig. 5 . Sensitivity of McCree equation plots [McCree, 1970, p. 227; our eqn 
(2)] to the principal parameters of the model (Fig. 1, Table lb). A 12-h day- 
length is applied followed by 12-h darkness. Light during the day was 
varied from near zero to 1000 J (PAR) nT 2 s 1 (4600 p-mol PAR nT 2 s -1 ). 
Plants are in 24 h exponential growth, attained after a few days from low 
initial values (Table la). The continuous line is for the default parameter 
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(C) senescence rate constant, k x ; (D) substrate utilization rate constant, k s ; 
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Fig. 6. The model is run in constant light flux densities (daylength = 24 h) 
from near zero to 1000 J (PAR) nT 2 s 1 (4600 (jimol PAR m~ 2 s _1 ) to a 
steady state with default parameter values (Table lb). Instantaneous values 
are plotted at daily intervals. (A) Leaf area index, L AI [eqns (3)] [black 
circles show inflexion points for total mass, eqns (18)]; (B) C substrate 
fraction, C s [eqn (11)]; (C) ratio of respiration to gross photosynthetic rate, 
r R:P [eqns (21)]. 



viewed as a 'maintenance' coefficient. According to the 
second of eqns (56), it depends on growth efficiency Y a , and 
is proportional to senescence rate k x and fraction of senescing 
material recycled f rec (Fig. 1). This differs from our extracted 
maintenance coefficient in eqns (31) and (42), which for 
high k s depends only on k x and Y G . 
If we re- write McCree's eqn (2) as 



Pnieht 



= ^'^n.day + d 



(57) 



using specific night-time respiration and specific daily net 
photosynthesis as defined in eqns (48) (these are measured 
quantities), then a replot of Fig. 5 using these changed vari- 
ables gives approximate parameters of 



f « (1 -/di)(l - Ya), J « fet(l -/di)(l - Y G )f K 



(58) 



Now fractional daylength [f dl , 
role [cf. eqns (56)-(58)]. 



eqn (43)] plays an important 



Steady-state simulations 

Short-term experiments (e.g. McCree, 1970) give valuable 
data, although interpretation may not be straightforward. 
However, the steady state of ecosystems is often of interest. 
Now we examine how the model parameters contribute to 
steady-state properties [eqns (36)-(42)]. First a true steady 
state when light flux density is constant is examined. 

Figure 6A shows how asymptotic L A1 is increased by light 
flux density. Inflexion points are shown; the ratio of inflexion 
L A1 to asymptotic L A1 decreases from 0-39 at low light (7 0 = 10 
W PAR m" 2 = 46 ujnol PAR m" 2 s" 1 ) to 0-16 at high light 
(/n= 1000 W PAR m" 2 = 4600 |xmol PAR m" 2 s" 1 ). At 
high light the growth curve becomes less 'logistic-like' 
(where the ratio is 0-5, e.g. Thornley and France, 2007, 
equation 5-19, p. 144). 

Figure 6B illustrates how C substrate fraction Cs [eqn (11)] 
has a value which, during exponential growth, depends on light 
flux density [eqns(27)], but which then decreases asymptoti- 
cally to a common value of 0-05333 [eqn (42)]. 

The ratio of respiration to gross photosynthesis [eqn (21)] is 
drawn in Fig. 6C. For young exponentially growing plants this 
is light dependent, with a low value at high light - decreasing 
light flux density increases this ratio [eqns (34)], until, at zero 
exponential growth rate, the ratio is 0-4 [eqn (35)]. In the 
steady state, a ratio of 0-4 [eqn (42)] is approached for all 
light flux densities even although the L A1 asymptotes are 
very different (Fig. 6A). This suggests that plant ecosystems 
not in the steady state can be expected to have differing 
ratios of respiration to gross photosynthesis, but in the steady 
state these tend to the same value [eqn (42)], which may be 
substantially independent of environment and ecosystem. 



12-hour daylength simulations 

If the system is driven by 12 h daylength (with daytime light 
levels doubled), the results are qualitatively as in Fig. 6 for a 
constant environment (now extracting and comparing daily 
averages or integrated fluxes). However, equilibrium L A1 
valuess are about half of the value in continuous light 
(Fig. 6A). Substrate fractions are much lower (about half) 
than in Fig. 6B during exponential growth. Surprisingly, the 
steady state daily average substrate fraction for the diurnal situ- 
ation is the same as the steady state value for the constant 
environment [Fig. 6B; eqns (42)]. The ratio of respiration to 
gross production as defined by eqn (46) ranges from 
approx. 0-25 (exponential growth, depending on light level) 
to 0-4 (steady state, for all light levels), showing the impor- 
tance of growth conditions and position in the life cycle 
between exponential growth and the steady state for this ratio. 



Respiration is zero during the daylight hours 

Next we explored the consequences of assuming that, during 
the 12 h light period of the diurnal environment, growth effi- 
ciency Y G is unity, instead of 0-75 [eqn (13)]. When illumi- 
nated, there is excess energy currency present allowing 
synthesis to occur without cost - there is no need to degrade 
photosynthetically fixed C to provide energy units for syn- 
thetic processes (Kromer, 1995). This assumption is made 
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for all light levels applied. These simulations are qualitatively 
as in Fig. 6. Compared with the 12 h day simulations reported 
in the last section, L A1 s are higher; substrate fractions are 
slightly lower (due to the greater dilution from increased 
growth - in fact there are higher amounts of substrate, M cs , 
Fig. 1). The ratio of 24 h respiration to daily gross photosyn- 
thesis [eqn (46)] is lower because daytime respiration is zero 
when Y G = 1 [eqns (13)]. For the respiration to gross photo- 
synthesis ratio, all curves approach the common asymptote 
°f '"R24h:/ , gday = 0-22 as compared with 0-4 in the last two sec- 
tions (e.g. Fig. 6C). This simulation demonstrates that any 
direct coupling of photosynthetically produced energy to 
growth processes could have large effects on respiration. 

Respiration in the dark 

Placing a plant, or part of a plant, in the dark, sometimes for 
prolonged periods, is sometimes used to give an estimate of 
'maintenance' respiration (Gary et ai, 2003). This topic has 
been partially addressed in the exponential growth section 
and eqns (28), (29) and (32). Here, in order to give some com- 
parison with Gary et ai (2003, fig. 1), the plant is grown in low 
and high continuous light [I 0 = 20, 1000 J PAR m" 2 s" 1 (92, 
4600 |xmol PAR m~ 2 s _1 ); Fig. 6] to a steady state, achieved 
by time t = 1000 d [eqns (42)]. Then the plant is placed in 
darkness and the time course of the specific respiration rate 
p, C substrate fraction Cs and plant mass M c is examined 
[eqns (14), (11), (18)]. Results are shown in Fig. 7. 

Plant mass is much higher at the higher light level, although 
the specific respiration rate p(ss) and substrate fraction Cs(ss) 
[eqns (42)] are independent of light level. After switching 
off the light (at time / = 1000 d), both these intensive quan- 
tities decline identically over some 10 d to their new constant 
values in the dark [eqns (32), (28)]. 

Comparison of Fig. 7 with fig. 1 of Gary et al. (2003), who 
plot tomato whole-plant respiration rate in prolonged darkness, 
shows that many extensions would be needed in order to simu- 
late real data. Issues include: multiple substrates, multiple 
structural components with different growth/senescence/ 
recycling characteristics, and possibly reaction switches in 
order to obtain increases in respiration rate. 



CONCLUSIONS 

The recycling model of Fig. 1 is able to account for three 
issues in plant respiration. 

First, required for the model's credibility, is to give an 
account of McCree's equation [McCree, 1970; our eqns (1) 
and (2)], which is measurement based, and arguably is why 
many researchers follow GMRP. This is spelt out in Fig. 4A, 
where it is assumed that plants are in exponential growth. 
The simulations of Fig. 5 reveal how McCree's parameters 
can be approximated by parameters of the model (Fig. 1) as 
given in eqns (56). Indeed our analysis suggests that the MR 
rate is not directly related to McCree's c parameter 
[eqn (2)], but is better regarded as the specific respiration 
rate at zero growth rate [eqns (31), (42)]. If McCree's equation 
is modified so that it relates more directly to measured quan- 
tities as in eqn (57), then the modified parameters can be 
approximated by eqns (58), where daylength is now important. 
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Fig. 7. Respiration decline in the dark. Plants are grown to a steady state 
[eqns (36)-(42)] at high and low continuous light [eqns (4)], 7 0 = 1000, 20 J 
PAR nT 2 s -1 (4600, 92 u,mol PAR m~ 2 s _1 ). At time t= 1000 d, light is 
switched off. Plant mass M c is given in eqns (18), specific respiration rate p 
[kg C respired (kg total C) 1 d '] by eqn (14), and substrate fraction Cs [kg sub- 
strate C (kg structural C) _I ] by eqn (11). 



The second issue concerns the nature of MR, which has long 
been a controversial concept (Wohl and James, 1942). Many 
may doubt the existence of a separate process which could 
be dubbed 'maintenance'. It is difficult to represent an 
MR process in crop and ecosystem models directly and 
effectively - such efforts easily give rise to unacceptable 
behaviour or ad hoc assumptions being introduced into the 
model. The approach of Fig. 1 provides a mechanistic alterna- 
tive, which it is straightforward to extend to deal with more 
detailed aspects of plants and crops. 

The third and last issue concerns the ratio of respiration R to 
gross photosynthesis P g [eqn (21)]. The topic has been 
addressed frequently over many years (e.g. Gifford, 2003) 
and recently by van Oijen et al. (2010). The model demon- 
strates that for a growing crop the ratio is expected to vary 
with time, as illustrated in Figs 2B and 6C. During exponential 
growth [eqns (34)], the ratio depends on light flux density, I 0 
[via eqns (8)]. In equilibrium [e.g. a mature plant ecosystem, 
eqns (42); see also eqn (35)] the ratio is constant, depending 
only on two parameters, and is independent of environment. 
This is demonstrated by the common asymptotes in Fig. 6C 
where the ratio is 0-4 for all values of light flux density. The 
value of 0-4 is likely to be conservative as growth yield 
1g = 0-75 is a well-established value for plant ecosystems 
where cellulose or similar compounds are a major component, 
and, as mentioned after eqns (16),/ rec = 0-5 may be a reason- 
able guess at the average fraction of senescing material 
recycled for real variable environments (note that / rec must 
lie between zero and unity). 

Gifford (2003) remarked that 'Plant respiratory regulation is 
too complex for a mechanistic representation in current terres- 
trial productivity models for C accounting...'. That may be 
true, but surely it is a council of despair - attempting to 
grasp and pin down complexity is often the first step to 
finding a way through the labyrinth. Many attempts have 
been made to understand in a simple manner the relationships 
between photosynthesis, growth, respiration and net pro- 
duction. These have not been successful, sometimes because 
they ignore some of the basic principles of reductionist 
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science and of explanatory dynamic simulation, promulgated 
by de Wit (1970) and others. 
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